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Abstract 

The capacitance between two adjacent nodes on an infinite square grid of 
identical capacitors can easily be found by superposition, and the solution 
is found by exploiting the symmetry of the grid. The mathematical 
problem presented in this work involves the solution of an infinite set of 
linear, inhomogeneous difference equations which are solved by the 
method of separation of variables. 
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1. Introduction 

The electric circuit networks has been studied well by Kirchhoff s 1 
more than 150years ago, and the electric-circuit theory is discussed in 
detailed by Van der Pol and Bremmer 2 where they derived the resistance 
between any two arbitrary lattice sites. Later on, at the ends of the 
eighties of the last century the problem is revived by Zemanian 3 , where 
he investigated the resistance between two arbitrary points in an infinite 
triangle and hexagonal lattice networks of identical resistor using Fourier 
series. For hexagonal networks, he discovered a new method of 
calculating the resistance directly from the networks. 

The problem is studied again by many authors 4 11 . For example, 
Cserti 8 and Cserti et. al 9 introduced in their papers how to calculate the 
resistance between arbitrary nodes for different lattices where they 
presented numerical results for their calculations. In recent works, we 
used Cserti 's method to calculate theoretically the resistance between 
arbitrary sites in an infinite square and Simple Cubic (SC) lattices and 
experimental comparison with the calculated values are presented 1011 . 
Finally, Wu 12 studied the resistance of a finite resistor network where the 
resistance between two arbitrary nodes is obtained in terms of the 
eigenvalues and eigenfunctions of the Laplacian matrix associated with 
the finite network. 

Little attention has been paid to infinite networks consisting of 
identical capacitances C. Van Enk 13 studied the behavior of the 
impedance of a standard ladder network of capacitors and inductors 
where he analyzed it as a function of the size of the network. In this paper 
we investigated analytically and numerically the capacitance between 
arbitrary lattice sites in an infinite square grid using the superposition 
principle. Also, the asymptotic behavior is studied for large separation 
between the two sites. An investigation of infinite complicated lattices 
and of lattices with missing capacitor (bond) is in progress. 

The physical situation is illustrated in Fig. 1. An infinite number of 
identical capacitors of capacitance C are connected to form an infinite 
square grid. The problem is to find the capacitance between arbitrarily 
spaced nodes. The basic approach used here is similar to that used by 
Paul 14 . 

Let a charge Q enter the grid at a node r and let it comes out of the 

grid at a distant point. Removing the return point to infinity then the 
problem is invariant under 90 "rotation, so the charge flowing through 
each of the four capacitors connected to the node will be equal. Therefore 
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each one of them will carry a charge of ^ . Thus, the resulting voltage 

drop between node r and an adjacent node r will be — . 

4 C 

Now, consider the case where a charge Q entering the grid at a distant 
point and exiting at the adjacent node, r . Again, the charge flowing will 

be — , and the voltage drop from r to r will be given by — . 
4 4C 

The superposition of the above two problems results in a new problem 

where a charge Q entering the node r and exiting the adjacent node r 

with a net voltage drop equal to — . The distant point is eliminated since 

the net charge there is zero, therefore the capacitance between adjacent 
nodes is2C. 

2. Node Voltages Analysis 

Consider an infinite square network consisting of identical 
capacitors such as that shown Fig. 1. Let the nodes be numbered from 
minus infinity to plus infinity in each direction, and let the voltage at the 
node (m,n) be denoted by V mn . Applying Kirchhoff s laws at node (m,n). 

Thus, one may write: 

Q m ,n = (V m ,n ~ V m ,n+\)C + (V mn — V m n _ l )C + (V mn — V m+l n )C + (V m n — V m _ l n )C . (1) 
If at node (m,n) the charge equal to zero, then Eq. (1) reduces to: 

Now, let a charge Q enter the node (0,0) and leave at infinity. Then 

Q mn = Unless m= and n= 0, Q 0fi =Q. (3) 
Or, we may write Eq. (3) as: 

f Q, m=n=0 

Q m , n =\ (4) 
I zero, otherwise. 

Equations (1 and 2) are the finite difference equivalent of Poisson's and 
Laplace's equations, respectively. 

3. Solution by Separation of Variables 
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Although the method of separation of variables is commonly 
applied to partial differential equations having suitable boundary 
conditions, it is equally applicable to difference equations 15 . 

Consider 

v m ,„ =QW(ma + in/3). (5) 
Substituting Eq. (5) into Eq. (2), we obtained the following: 

4 exp(ma + in/3) = exp((m + \)a + in/3) + exp((w - \)a + in/3) + 

exp(ma + i(n + l)/3) + exp(ma + i(n - 1)/?). (6) 

The above equation can be rewritten as: 

4 exp(ma + in/3) = exp(ma + in/3)(exp(a) + exp(-a) + exp(i/3) + exp(-//?)) . (7) 
Thus, Eq. (2) is satisfied provided that: 

Cosha + Cos/3 = 2. (8) 

Noting that our aim is to solve the problem with a source at (0,0), this 
implies that: 

v -V =v -V -V - (9) 

m,n n,m —mji m,-n m,-n V / 

Therefore, take 

V mn - exp(-\m\a)Cosn/3 + exp(-\n\a)Cosm/3 . (10) 

The above functions V mn do not satisfy the source- free difference 

relation given by Eq. (2) along the lines n=0 and m=0. This is due to the 
absolute value sign in the exponential terms, so there will be residual 
charges entering or leaving the grid at each node along these lines. 

To find the external charges Q mn which produce the voltage 

pattern V m n , we may write from Eq. (1): 

= Ay _y _y _y _y ■ 

^ r o,o y o,l y o,-l y l,o r -\,oi 

= 8 - 4(Cos/3 + exp(-a)) . (11) 
Using Eq. (8), then Eq. (11) can be simplified as: 
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Q^ = 4Sinha. (12) 
C 

In a similar way, for n * 



C 



= 4V — V -V -V -V 

^ r o,n v o,n+l v o,n-l V 1,11 r -1,11 1 



= ICosnpil - Cos/3 - Cosha) + 2 exp(\n\a) {2 - Cos (5 - Cosha} . (13) 
Again, using Eq. (8) we can simplify Eq. (13) as: 

Q^ = 2Cosn/3Sinha. (14) 

4. Charge Entering at (0,0) 

The assumed node voltage pattern V m n (/?) requires external charges 

Q m „ not only at node (0,0) but at all nodes for which either m= or n= 0. 

Thus, it is necessary to form a superposition of such voltages with 
different values of/?, to suppress all external charges except the one at 

(o,o;. 

Now, let V m n (/?) having the following form: 

K, n =] F(PW m , n {P)dp. (15) 

o 

where the limits of the integral have been chosen to cover the entire 
applicable range of values of p. The function F(/?)is an amplitude 
function that must be chosen to make Q = Q and Q o n = , when n * . 
Substituting Eq. (15) into Eq. (1), we obtained: 

^ = ]F(P)Q m ,AP¥P. (16) 

^ 

From Eqs. (12 and 13), we have: 

O 71 

^ = f [F(P)ASinhaW . (17) 

C o 

and 

= \[F(P)2SinhaCosnpyip . (18) 
C { 
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where Q o n = Q a _ n = Q n o = Q_ n o .(i.e. Q o n is symmetric). 
The expression for F(j3) can be obtained by inspection as: 

F(P) = & . (19) 

C(4nSinha) V 7 

which satisfies the above condition. Thus, Eq. (15) becomes: 







Q r (exp(- \m \a)CosnP + exp(- « \)CosmJ3) 

— — dp . (20) 

AnC { Sinha 



o 



Note, that a is a function of p . They are related by Eq. (8). 
5. Capacitance between Two Points in a Large Grid 

We mentioned earlier that the capacitance between (0,0) and (m,n) 
could be obtained directly from the solution of the problem in which the 
charge Q enters at (0,0) and leaves at infinity. In terms of the node 
voltages, the capacitance C m n can be written as: 

C mn = • (21) 

2(V -V ) 

V o,o m,n ' 

Using Eq. (20), one may write Eq. (21) as: 



Cm ' n H 2 " ex P(- 



m 



2kC 

a)CosnP — exp(- n a)CosmP) 
Sinha 



(22) 



dp 



Or, one may write C m n (see Appendix A) as: 

C mn = — . (23) 

| (1 - exp(- n a)Cosm P) 

{ Sinha 



The integrals in Eqs. (22 and 23) have to be evaluated numerically. 
Results of C m n for values of (m,n) ranging from (0,0) - (10,10) are 

presented in Table 1 below, and here are some special cases: 

(a) C 0;0 
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Using eq. (22) with n=m=0. Then 

r (2 - exp(-0)CosO - exp(-O)CosmO) 



C = — = = oo . (24) 

"" (2 exp( ())C n 

Sinha 



as expected. This can be explained as a parallel capacitance with zero 
separation between its plates. 

(b) c 04 

Again, using Eq. (22) with n=0, m=l. Thus 

C o j = — = 2C . (25) 

| (2 - Cos/? - exp(-oQ) 

* Sinha 
where we have used 2 - Cos/3 - exp(-a) = Sinha . 

The above result is the same as mentioned at the end of the introduction. 
We may explain this result as two identical capacitors connected together 
in parallel. 

(c) Asymptotic form for large morn 

When either m or n are large, the exponential terms given in Eq. (23) 
become negligible except when a is very small. When a is very small; 

a 2 

Cosj3 = 2-Cosha*\- — . (26) 
So that 

a « Sinha - J3 . (27) 
Suppose m is large, then Eq. (23) can be rewritten as: 

= 7V[ ; ; 1 h 



C U\-exp(-mj3)Cosnj3) U\/Siha -\) 



1 1 -]• (28a) 



? exp(- m /3)Cosn/3 ? exp(- m /3)Cosn/3 



^(f + f + f)- (28b) 

^2 ^2 ^3 
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The first integral can be expressed in terms of the exponential integral 
Ein(z) 16 



= |(l-exp(-A/^))^ (29) 

P 

So, /, = Re{Ein[ft(n + im)]} . 

The second integral can be integrated numerically, and it is found to be 
1 2 = -0. 1049545 . In the third integral the exponentials are negligible except 
for small values of a and/? , and for those values a « Sinha « p . So, 7 3 can 
be neglected. Thus, Eq. (28) becomes: 

C ~ Re{Eiri7r(n + im)]} 0.1049545' 

For large values of its argument, Eiriz)^- In z + 0.57721. Therefore, Eq. 
(30) can be rewritten as: 

m,n /T 1 \ 

C ~21n(n 2 +m 2 ) + 4\n ^ + 0.22038855' 

For reasonable values of m and n, the asymptotic form (i.e.: Eq. (31)) 
gives an excellent approximation, and from this equation we can show 

C 

that as any of m and n goes to infinity then -^--»0, which can be 

explained as a parallel capacitance with infinite separation between its 
plates. Finally, it is clear from Eq. (23) that C mn = C_ m _ n which is 

expected due to the inversion symmetry of the infinite square grid. 

6. Results and Discussion 

In this work, the capacitance between the site (0,0) and the site 
(m,n) ; in an infinite square grid consisting of identical capacitors is 
calculated using the superposition of charge distribution. The capacitance 
C m n is expressed in an integral form which can be evaluated numerically 

or analytically. 
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The asymptotic form for the capacitance as m or/and n goes to 
infinity is investigated where it is shown that it goes to zero. 

In Figs. (2-5) the capacitance is plotted against the site (m,n). Fig. 
(2) shows a three dimensional plot of the capacitance as a function of m 
and n. One can see from the figure that as m or/and n increases then C m n 

decreases up to zero at infinity as expected before (i.e. see Eq. (31)). 

Figures (3-5) show the capacitance C m n as a function of the site 

(m,n) along the directions [10], [01] and [11]. From these figures we can 
see that the capacitance is symmetric along these directions, and this is 
due to the inversion symmetry of the infinite square grid. Also, the 
figures show how the capacitance C m n goes to zero as any of m or n goes 

to infinity. 
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Figure Captions 



Fig. 1 An infinite number of identical capacitors of capacitance C are 
connected to form an infinite square grid. 

Fig. 2 The capacitance C m n in terms of m and n for an infinite square grid. 
Fig. 3 The capacitance C m n in terms of the site along [10] direction. 
Fig. 4 The capacitance C m n in terms of the site along [01] direction. 
Fig. 5 The capacitance C m n in terms of the site along [11] direction. 
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C(m,n)/C 



Fig. 2 
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Table Captions 

Table 1: Numerical values of C m n in units of C for an infinite square grid. 



Table 1 









( m 


L m,» /C 




. 822545 


(9,1) 


0.907753 


(6,2) 


oo 


(0,0) 


0.818628 


(9,2) 


0.892285 


(6,3) 


2 


(1,0) 


0.812497 


(9,3) 


0.874193 


(6,4) 


1.5708 


(i,i) 


0.804631 


(9,4) 


. 85517 


(6,5) 


1.37597 


(2,0) 


0.795539 


(9,5) 


0.836326 


(6,6) 


1 .29326 


(2,1) 


.785687 


(9,6) 


0.882207 


(7,0) 


1.1781 


(2,2) 


0.775459 


(9,7) 


0.879628 


(7,1) 


1.16203 


(3,0) 


0.765148 


(9,8) 


0.872324 


(7,2) 


1.13539 


(3,1) 


.754964 


(9,9) 


0.861357 


(7,3) 


1 .08177 


(3,2) 


0.801699 


(10,0) 


0.847985 


(7,4) 


1.02443 


(3,3) 


0.800666 


(10,1) 


0.833344 


(7,5) 


1.04823 


(4,0) 


0.797649 


(10,2) 


0.818295 


(7,6) 


1.03649 


(4,1) 


0.792868 


(10,3) 


0.803421 


(7,7) 


1.00814 


(4,2) 


0.786636 


(10,4) 


0.850222 


(8,0) 


0.972869 


(4,3) 


.779303 


(10,5) 


0.848397 


(8,1) 


0.937123 


(4,4) 


0.771206 


(10,6) 


0.843152 


(8,2) 


0.974844 


(5,0) 


0.762645 


(10,7) 


0.835079 


(8,3) 


0.968523 


(5,1) 


0.753862 


(10,8) 


0.824942 


(8,4) 


0.951831 


(5,2) 


.745047 


(10,9) 


0.813496 


(8,5) 


0.929041 


(5,3) 


0.736338 


(10,10) 


0.801381 


(8,6) 


. 90391 


(5,4) 






0.789079 


(8,7) 


. 878865 


(5,5) 






0.776929 


(8,8) 


0.922313 


(6,0) 






0.823894 


(9,0) 


0.918443 


(6,1) 
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Appendix A 

Instead of using the functions V m n defined in Eq. (10), let us use the 
following functions 

W mn = exp(-\m\a)Cosfi . Al 

Equation (Al) is a source free everywhere except a long the line m=0. 
Thus, the external charge Q 0n can be written as: 

= 2Cosn /3Sin h a . A2 
Now, let V m n be given as: 

V m , n =\F(P)W m , n (P)dp. A3 

o 

The corresponding expression for the external charge is 

— ^_ = J F(j3)2SinhaCosn J3dj3 . A4 
C 

Using Fourier cosines series, one can write (using Eq. (A4)) 

00,0 +^TjQ^ C0Sn P 

27iF(j3)Sinha = = . A5 

For our case considered here 

2nF{P)Sinha = -. A6 
Or, we may write 

F(J3) = ^ . A7 

2nCSinha 

Substituting Eqs. (Al and A7) into Eq. (A3), we obtained: 
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O r exr>(-\m\a)Cosn/3 

v m n = g r FV 1 1 ' —dp . 

2nC * Sinha 

Finally, the capacitance C m „ can be obtained by inserting Eq. (A8) into 
Eq. (21). Thus we get: 



A8 



C. 



c 



r 1 - exp(- m a)CosnB 

| \ -<//? 



A9 



Sinha 
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